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We propose an extension of the hard- wall AdS/QCD model by including the Chern-Simons term 
required to reproduce the chiral anomaly of QCD. In the framework of this holographic model, 
we study the vertex function F 1T1 * 1 * (Qi, Ql) which accumulates information about the coupling 
of the pion to two (in general virtual) photons. We calculate the slope of the form factor with 
one real and one slightly virtual photon and show that it is close to experimental findings. We 
analyze the formal limit of large virtualities and establish that predictions of the holographic model 
analytically (including nontrivial dependence on the ratio of photon virtualities) coincide with those 
of perturbative QCD calculated for the asymptotic form of the pion distribution amplitude. We also 



investigate the generalized VMD structure of 
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■ (Qi, Ql) in the extended AdS/QCD model. 
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I. INTRODUCTION 



The form factor F 7 » 7 » w o (Q 2 , Q\) describing the cou- 
pling of two (in general, virtual) photons with the light- 
est hadron, the pion, plays a special role in the studies of 
exclusive processes in quantum chromodynamics (QCD). 
When both photons are real, the form factor F 7 » 7 . 7r o (0, 0) 
determines the rate of the 7r° — ► 77 decay, and its value 
at this point is deeply related to the axial anomaly 
Because of this relation, the 7*7*71-° form factor was an 
object of intensive studies since the 60's @-@|- 

At large photon virtualities, its behavior was stud- 
ied @, H, H within perturbative QCD (pQCD) fac- 
torization approach for exclusive processes 0, [ToL [ll[. 
Since only one hadron is involved, the 7*7*71"° form fac- 
tor has the simplest structure for pQCD analysis, with 
the nonperturbative information about the pion accu- 
mulated in the pion distribution amplitude ^^(x) intro- 
duced in Refs. [12L fl3 ]. Another simplification is that 
the short-distance amplitude for 7*7*71-° vertex is given, 
at the leading order, just by a single quark propagator. 
Theoretically, most clean situation is when both pho- 
ton virtualities are large, but the experimental study of 
F 7 . 7 . 7r o(Q 2 , Ql) in this regime through the 7*7* — * 7r° 
process is very difficult due to very small cross section. 

The leading-twist pQCD factorization, however, works 
even if one of the photons is real or almost real. Further- 
more, this kinematics is amenable to experimental inves- 
tigation through the 77* — > 7r° process at e + e _ colliders. 
Comparison of the data obtained by CELLO [l4[ and 
CLEO [ll| collaborations with the original leading 0, [1] 
and next-to-leading order 

PS EMM! pQCD predic- 
tions amended by later studies[20l l2ll| that incorporate 
a more thorough treatment of the real photon channel 
using the light-cone QCD sum rule ideas provided im- 
portant information about the shape of the pion distri- 
bution amplitude tp v (x) (for the most recent review see 



|22j|). The momentum dependence of F 7 » 7 »„.o(0, Q 2 ) form 
factor was also studied in various models of nonpertur- 
bative QCD dynamics (see, e.g., [23j]-[37| and references 
therein). 

In the present paper, our goal is to extend the holo- 
graphic dual model of QCD to incorporate the anomalous 
F 7 * 7 » T o (Qi, Q 2 .) form factor. During the last few years 
applications of gauge/gravity duality [38| to hadronic 
physics attracted a lot of attention, and various holo- 
graphic dual models of QCD were proposed in the lit- 
erature (see, e.g., [3!|-[5!|)- These models were able to 
incorporate such essential properties of QCD as confine- 
ment and dynamical chiral symmetry breaking, and also 
to reproduce many of the static hadronic observables (de- 
cay constants, masses), with values rather close to the 
experimental ones. 

As the basis for our extension, we follow the holo- 
graphic approach of Refs. [43|, [45|, and then intend to 
proceed along the lines of formalism outlined in our re- 
cent papers [6(| [fiTIi where it was first applied to form 
factors and wave functions of vector mesons [63, HH (ten- 
sor form factors of vector mesons were considered in [63]) 
and later [6l[ to the pion electromagnetic form factor. 
However, a straightforward application of the approach 
of Refs. [H EE IH !l| to F^.^iQhQl) form factor 
gives a vanishing result. There are two obvious reasons 
for such an outcome. 

First, the five-dimensional (5D) gauge fields B a (x, z) in 
the AdS/QCD lagrangian of Refs. [U, [4{| are only dual 
to the 4D isovector currents J a (x). On the other hand, 
a nonzero result for the matrix element (0| Jg M </emK ) 
defining the i ? 7 * 7 * 7r o(Qi, Ql) form factor may be ob- 
tained only when the electromagnetic currents J^ M , Jem 
have both isovector and isoscalar components, which is 
the case in two-flavor QCD, but not in the holographic 
models of Refs. (HHl]. 

Thus, we need an AdS/QCD model that includes gauge 
fields in the 5D bulk which are dual to both isovector and 
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isoscalar currents. The natural way to do this is to extend 
the gauge group SU(2) L ® SU(2) R in the bulk up to 
U(2)l®U(2)r. After explicit separation of the isosinglet 
and isovector parts, the new 5D field can be written as 



B^ = t a B° 



t • 



■j-, with 1 being the unity matrix. The 

B part is dual to 4D isosinglet vector current. 

Even after this modification, the AdS/QCD action 

gives zero result for the correlator (0|4 /_1} J„ {/ ~ 0} J^|0) 

and the ax- 



involving two vector currents 



{/=!} 



.{1=0} 



ial current (that has nonzero projection onto the pion 
state). To bring in the anomalous amplitude into the 
model, the next step (similar to [67]]) is to add a Chern- 
Simons term [68| to the action. After these extensions, 
the calculation of the F 7 . 7 . w o (Ql, Q%) form factor may be 
performed using the methods developed in Refs. [60Ll6lj|. 

The paper is organized in the following way. We start 
by recalling, in Section II, the basics of the hard-wall 
model, in particular, the form of the action given in 
Ref. [HI and our results [^H for the pion wave func- 
tion. In Section III, we consider the generalization of 
the AdS/QCD model that includes isoscalar fields and 
Chcrn-Simons term. Using this extended model we de- 
scribe the calculation of the F 1 * J * 7T o(Ql, Q 2 ) form factor 
and express it in terms of the pion wave function and two 
bulk-to-boundary propagators for the vector currents de- 
scribing EM sources. In Section IV, we study the results 
obtained within the extended AdS/QCD model with one 
real and one slightly virtual photon and calculate the 
value of the Q 2 -slope of the form factor. We also discuss 



the formal limit of large photon virtualities, and compare 
these results to those obtained in pQCD. In Section V, we 
study the generalized VMD structure of the AdS/QCD 
model expression for the i 7 ' 7 » 7 * 7r o (Q\ , Q 2 ) form factor. Fi- 
nally, we summarize the paper. 



II. BASICS OF HARD- WALL MODEL 

In the holographic model, QCD resonances correspond 
to Kaluza-Klein type excitations of the sliced AdS5 space 
with the metric 



MjN 



gMNdx dx 



1 



{r\p, v dx^ dx v — dz 



(1) 



where rj^v = Diag (1, — 1, — 1, — 1) and [i. v = (0,1,2,3), 
M,N = (0, 1,2,3, z). The basic prescription is that 
there is a correspondence between the 4D vector and 
axial- vector currents and 5D gauge fields V^(x,z) and 
A a ^(x, z). Furthermore, since the gauge invariance corre- 
sponding to the axial- vector current is spontaneously bro- 
ken in the 5D background, the longitudinal component 
of the axial-vector field becomes physical and related to 
the pion field. 

A. AdS/QCD action 



The action of the holographic model of Ref. [4c 
be written in the form 



can 



J 



•^AdS — Tr 



d 4 a 



dz 



-(D M X)\D M X) + ^X 



*9P 



2 {Bf* T ?B 



(L) B (L)MN + 



r>MN R \ 
(R) ts (R)MN) 



(2) 



where DX = dX — iB (L) X 



iXB (R) , (B {LjR) =V±A) 



and X(x,z) — v(z)U(x, z)/2 is taken as a product of 
the chiral field U(x,z) = exp {2it a ir a (x, z)\ (as usual, 
t° = (T a /2, with a a being Pauli matrices) and the func- 
tion v(z) = (m q z + az 3 ) containing the chiral symmetry 
breaking parameters m q and er, with m q playing the role 
of the quark mass and a that of the quark condensate. 

In general, one can write A = A± +Ah, where A± and 
Ah are transverse and longitudinal components of the 
axial-vector field. The spontaneous symmetry breaking 
causes Am to be physical and associated with the Gold- 
stone boson, pion in this case. The longitudinal compo- 
nent is written in the form: 



A\ M {x,z) = d M tp a (x,z) 



(3) 



Then ip a (x,z) corresponds to the pion field. This Higgs- 
like mechanism breaks the axial- vector gauge invariance 
by bringing a z-dependent mass term in the A-part of 
the lagrangian. 



Varying the action with respect to the transverse gauge 
fields and A° gives equations of motion for these 
fields describing (via the holographic correspondence) 
the physics of vector and axial- vector mesons. Vari- 
ation with respect to A?, and A a z gives two coupled 
equations for the chiral field 7r°(x, z) and the pion field 
tjj a (x, z). It is convenient to work in Fourier representa- 
tion, where V^(p, z) — V^(p)V(p, z) is the Fourier trans- 
form of V^(x,z), and A^(p, z) is the Fourier transform 
ofA-(x,z). 



B. Vector channel 

The vector bulk-to-boundary propagator 



V(p, z) = g 5 



^ -pi + Ml 

m=l " l 



(4) 
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has the bound-state poles for p 2 = M 2 , with the res- 
onance masses M n = 7o,n/ z o (7o,n is n * h zer0 °f the 
Bessel function Jo (a;)) determined by the eigenvalues of 
the equation of motion 



An essential ingredient of form factor formulas is the 
vector bulk-to-boundary propagator J{Q,z) = V{iQ,z) 
taken at a spacelike momentum p with p 2 = —Q 2 . It can 
be written in a closed form as 



d z 



z 



+ l M 2 ^(z) = 0, (5) 



subject to boundary conditions V^(0) = d z ipX( z o) = 0. 
The eigenf unctions of this equation give the u ip n wave 
functions 



V2 



zJl{M n z) 



(6) 



for the relevant resonances. The coupling constants f n 
are determined from the "0 wave functions through 



fn 



1 

95 



\/2M„ 



In Ref. [60(| , we introduced 



z=0 95Z Jl(lO,n) 

<t> wave functions" 



(7) 



M„z 



d^X(z) 



V2 



ZoJl(lO,n) 



J (M n z) , (8) 



which give the couplings g^fn/Mn as their values at the 
origin, just like the [L = 0) bound state wave functions 
in quantum mechanics. Moreover, these functions sat- 
isfy Dirichlet b.c. <p!^{zo) = 0. Physically, the ip wave 
functions describe the vector bound states in terms of the 
vector potential V^, while the (f> v wave functions describe 
them in terms of the field-strength tensor d z V^ = V zfJi (we 
work in the axial gauge, where V z =0). 



J(Q,z) = Qz 



K 1 {Qz)+h{Qz 



K (Qz Q ) 
Ia{Qz ) 



(9) 



The function J{Q, z) satisfies the relations J(Q, 0) = 1, 
J(0, z) = 1 and d z J(Q, z ) = 0. 



Pion channel 



An important achievement of the hard-wall model of 
Ref. [43| is its compliance with the Gell-Mann-Oakes- 



Renner relation 



m q that produces a massless pion 



in the m q — limit. The fits of Ref. [43j give very small 
value m q ~ 2 MeV for the "quark mass" parameter m q , 
so it makes sense to resort to the chiral m q = limit, 
which has an additional advantage that solutions of equa- 
tions of motion in this case can be found analytically. 
The pion wave function ip( z ) is introduced through the 
longitudinal part of the axial-vector field: 



(10) 



The bulk-to-boundary part ir(z) of the 7r°(z) field, in the 
zeroth order of the of m 2 expansion proposed in Ref [43j 
tends to —1. Then the equation for ^(z) = ip(z) — ir(z) 
is exactly solvable, with the result 



*(z) = zr(2/3)(|) 



1/3 



I-1/3 {az 3 ) - h/3 



az 



h/3 



1-2/3 (otzl) 



(11) 



where a = g^cr/S (one may use here Airy functions in- axial- vector current 
stead of I-i/3(x),I 2 /3{x), cf. [HI)- The pion wave 



function ^(z) coincides with the axial- vector bulk-to- 
boundary propagator -4(0, z) taken at p 2 = 0. It is 
normalized by $(0) = 1, satisfies Neumann boundary 
condition ^>'(zq) = at the IR boundary and, due to the 
holographic correspondence, has the property that 



(12) 



For the neutral pion, it is convenient to define f T 
through the matrix element of the 173 projection of the 



J* V(p)> EE 



7r °b) ) = *JWv ■ 
(13) 



Then, the (experimental) numerical value of /„■ is 
92.4 MeV. Matching AdS 5 result £ AdS (p 2 ) - lnp 2 /(2 ff 2 ) 
and QCD result ^ CD (p 2 ) ~ lnp 2 /(87r 2 ) for the 
large-p 2 behavior of the correlator of jfi' 3 currents gives 
95 = 2tt [43j. Analyzing pion EM form factor, one deals 
with charged pions, and the choice of the axial-vector 
current as J^f = dj^su is more natural. Then 
fi c) = 130.7 MeV, while X^ D (p 2 ) ~ ^ P 2 /(4tt 2 ), hence, 
3g = \f2-K (6{|. Of course, the combination g 2 f 2 , be- 



4 



ing the ratio of the coefficient of the pion pole contribu- 
tion T. n (p 2 ) ~ fl/p 2 to the coefficient ~ 1/gf in £(p 2 )'s 
large-p 2 behavior, remains intact: 

g\fl = 4** ft = 2^(/M) 2 = (gff^ ^ s /2 , (14) 

where so f=s 0.67 GeV 2 . It is this convention- independent 
combination that enters Eq. (JT2J) . 

Again, it is convenient to introduce the conjugate wave 
function 1611: 



*(*) 



so \z 



(15) 



It vanishes at the IR boundary z = zq, i.e., 3>(#o) = 
and, according to Eq. (p~2|) . is normalized as $(0) = 1 at 
the origin. From $(0) = 1, it follows that the pion decay 
constant can be written as a function 



glfl = 3 . 2 1 / 3 r(2/3) 72/3 ^ 



v 2 /3 



T(l/3) I_ 2/3 (az 3 ) 



(16) 



of the condensate parameter a and the confinement ra- 
dius zo- Note that the magnitude of a is independent of 
the 55-convention, while the value of a depends on the 
35-convention used. 



After fixing zq 



zq - 
for 



= z? = (323 MeV) 
a = (424 MeV) 3 . 



through the p- meson mass, 
, experimental f v is obtained 
For these values, the argument 



ccZq of the modified Bessel functions in Eq. ([T 



equals 2.26 
then 



a . Since h/3{a)/I-2/ 3 (a) « 1 for o > 1, 



glfl « 3 • 2 1 / 3 



£(2/3) 

r(i/3) 



a 



2/3 



(17) 



i.e., the value of f„ is basically determined by a alone 
(th e same observation was made in the pioneering paper 
[45| and in a recent paper [6~i| in which the pion channel 
was studied numerically). 



In Fig.[TJ we illustrate the behavior of the pion wave functions \f r and <E> representing them W — » i/j(£, a), $ — > a) 
as functions of dimcnsionlcss variables C = z / z o an d a = az^. For a = (i.e. when the chiral symmetry breaking 
parameter a vanishes), the limiting forms are ip(£,0) = 1 and ip(C,0) = 1 — C 4 [61]. As a increases, both functions 
become more and more narrow, with ?/>(£, a) becoming smaller and smaller at the IR boundary £ = 1. 
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FIG. 1: Functions V(C,a 
a = 10 (lowermost lines). 



(left) and tp(£,a) (right) for several values of a: a — (uppermost lines), a = 1, a — 2.26, a — 5, 



III. ANOMALOUS AMPLITUDE 
A. Isosinglet fields 

The 7r°7*7* form factor is defined by 

J (7T,p|T{j£ M (x) J£ M (0)} \0)e-^d*x 

= e^q la q 2p F 1 , r ^{QlQl). i (18) 



where p = q\ + q 2 and q\ 2 = — Q 2 2 - Its value for real 
photons 

N 

F 7«7^o(0,0) = I ^- (19) 

is related in QCD to the axial anomaly. 

Since 7T° meson is described by the third component A 3 
of an isovector field, only the isovector component of the 
product of two electromagnetic currents gives nonzero 
contribution to the matrix element in Eq. (|18[) . In QCD, 
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the electromagnetic current is given by the sum 



tEM 



_ t{/=1},3 , 1 t{/=0} 

- j p + 3 J p 



(20) 



of isovector jff ^ ,a {x) ("p-type") and isosinglet 
Jp 1-0 ^ { x ) ("w-type") currents, 

^ /=1} ' 3 = ~ («7m« - htA = W^q . (21) 

4 I=0} = \ («7m« + <% rf ) = g^/* 1 * • 

As a result, the matrix element (tt°| J em J em |0) is 
nonzero since it contains (7T°| j{ /=1 }< 3 j{ /= °}|0) ~ 
Tr(i 3 i 3 ) parts. 

To extract JL* 7 » 7r o (Qf , Q2) form factor via holographic 
correspondence, we consider the correlator of the axial- 
vector current 



(22) 



and vector currents jj 1 I ~ 1 ^' b , J^ 7-0 ^. To proceed, we 
need to have isoscalar fields on the AdS side of the holo- 
graphic correspondence. This is achieved by gauging 
U(2) L <g> U(2) R rather than SU{2) L ® SU(2) R group in 
the AdS/QCD action, i.e., by substituting i a B£ by 



(23) 



in Eq. ([5]). Then the 4D currents correspond to the fol- 
lowing 5D gauge fields 

J$< a (x)^A;(x,z) , (24) 

jy =i} ' a (x)^v^ x ,z) . 



B. Chern-Simons term 

The terms contained in the original AdS/QCD action 
([2|) cannot produce a 3-point function accompanied by 
the Levi-Civita e^ af3 factor. However, such a contribu- 
tion may be obtained by adding the Chern-Simons (CS) 
term. We follow Refs. [65|, [|56| in choosing the form of 
the 0(B 3 ) part of the D = 5 CS action (the only one we 
need). In the axial gauge B z = 0, it is written as 

4s[£] = k ^ e ^ PaTl J d 4 xdz(d z Bp) 



B v T t 



pa 



(25) 



where T is the field-strength generated by B and k should 
be an integer number adjusted to reproduce the QCD 
anomaly result. Then, in the U (2)j, (8 U(2)r model, the 
relevant part of the CS term is [67[ 



Taking into account that Bl.r — V ± A, and keeping 
only the longitudinal component of the axial- vector field 
A — > A\\ (that brings in the pion), for which F^ v = 0, we 
have 



^ S ds = fc^Wd** r°d Z 



(27) 



(d p VZ) A\ a d z V v ) + (d p %) At d z V v a 



where d z = d z — d z . Representing A^ a = d a ip a and inte- 
grating by parts gives 



oAdS 
°CS 



k^ Vpa ld 4 xl dz 



(28) 



i{d z r) {d p V«) (d a V v )-il> a d z {d p V^d a V v 



Note that the z-derivative of -0° is proportional to the 
function of the pion (|15p . which, as we argued in 
Refs. [iil [ML H3|, is the most direct analog of the 
quantum- mechanical bound state wave functions, i.e., it 
is the derivative d z ip a (z) rather than ip a (z) itself that is 
an analog of the pion 4D field. Then, the first term in 
the square brackets in Eq. (f2"5)l has the structure similar 
to the TTLup interaction term 



C^ p = ^f-e^ir" («) {d a oj p ) , (29) 



obtained in the hidden local symmetries approach [701 ] 
from the anomalous Wess-Zumino [7l|,[72j lagrangian (see 
also a review [73j ) . Here g is the universal gauge coupling 
constant of that approach. 

Integrating by parts over z, the second term of Eq. (j2"5|) 
can be also converted into a contribution of this form plus 
a z = zq surface term: 



cAdS 



24ir 2 

20 



dz (d z ^ a ) (d p vz) (dX) 



(30) 



The latter can be eliminated by adding a compensating 
surface term into the original Chern-Simons action, so 
that the resulting anomalous part of the action in the 
extended AdS/QCD model 



^anom = fc o ^upa f fa 

8n 2 



dz 



{d z r) (d P V«) (d„%) (31) 



Scs b l B L, B R ] - S^ s [B L ] ~ S'cs • ( 26 ) has the structure of Eq. flU}. 



?(3) 
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C. Three-point function 



The action (|30|) produces the 3-point function 

<^< 3 (-^H^ EM te)) 

= T a ^ v (p,q 1 ,q 2 ) i(27r) 4 <5 (4) ((7i + q 2 - p) 



with 



12tt z p z 

Here, p is the momentum of the pion and 91,92 are the 
momenta of the photons. 

The "bare" function K b (Q\, Q 2 ) is given by 



K b (QlQ 2 



k f z ° 

2; 77 / J{Qi,z)J(Q 2 , z)d z *P(z) dz ,(33) 



2 7o 



where J(Qi,z), J{Q 2 ,z) are the vector bulk-to- 
boundary propagators (O, and ip(z) is the pion wave 
function (fTTj) . 



JziQi = 0,zo) = at the real photon point. On the 
other hand, S(Qi = 0, zq) — 1, and our final model for 
F^ 7 « 7 . (Qf , Q 2 ) corresponds to the function 



K{Q\,Qt) = ^(z )J(Q 1 ,z )J(Q 2 ,z ) 



(36) 



J(Q 1) z)J{Q 2 ,z)d z *(z)dz 



The extra term provides K(Q,0) — 1, and since 
J(Q,z Q ) ' 



Io(Qzo) 



(37) 



it rapidly decreases with the growth of Qi and/or Q 2 . 
We will see later that in these regions the behavior of 
F ff7 . 7 . (Qf, Q 2 ) is determined by small-z region of inte- 
gration. Thus, the effects of fixing the ^(zq) ^ artifact 
at the infrared boundary are wiped out in the "short- 
distance" regime. 



IV. MOMENTUM DEPENDENCE 



D. Conforming to QCD axial anomaly 



Small virtualities 



For real photons, i.e., when Q\ = Q\ = 0, the value 
of the F 7r7 * 7 » (Qi, Q 2 ) form factor in QCD (with mass- 
less quarks) is settled by the axial anomaly, which cor- 
responds to ifQ CD (0, 0) = 1. Our goal is to build an 
AdS/QCD model for F 7r7 » 7 . (Q\, Q^) that reproduces this 
value. Taking Q\ = Q\ — 0, we have J{Q,z) = 1, and 
Eq. {33]) gives 



k f z ° 

K b (0,0) = -- d 7 Mz)dz 
* Jo 



2 V>Oo) 
1 - tt(z ) 



(34) 



On the IR boundary z 
from Eq. (fTTj) is 



Zq 



the pion wave function \I/(z) 



V3r ( 2/3) / J_ 
( 0) 7r/_ 2/3 (a) U 2 



1/3 



(35) 



As we discussed above, experimental values of m p and 
f„ correspond to a — 2.26, which gives ^(20) = 0.14. 
The magnitude of ^S(zq) rapidly decreases for larger a 
(e.g., ^{z )\ a=4 « 0.02, see Fig. [JJ). Still, the value of 
^(zq) is nonzero at any finite value of a, and it is im- 
possible to exactly reproduce the anomaly result by sim- 
ply adjusting the integer number k. To conform to the 
QCD anomaly value K^ CT> (0, 0) = 1, we add a surface 
term compensating the ^(zq) contribution in Eq. (34) 
and then take k = 2. To fix the form of the surface 
term, we note first that it should have the structure of a 
VVA 3-point function taken on the z = zq surface. Fur- 
thermore, using the derivatives 3' z {Qi, zq) of the bulk-to- 
boundary propagators in this term is excluded, because 



If one of the photons is real Q\ = 0, while another is 
almost real, Q\ — Q 2 -C l/zfi, we may use the expansion 



J(Q,z) 



\Q 2 z 2 



In 

4 



0(Q 4 



(38) 



which gives (for zq = Zq and a — clq = 2.26) 



K(0,Q 2 ) ~1 - 0.66 



Q 2 4 



4 

)2 



1 - 0.96- 



(39) 



The predicted slope is very close to the value 1/m 2 ex- 
pected from a naive vector-meson dominance. Experi- 
mentally, the slope of the F 7r7 * 7 * (0, Q 2 ) form factor for 
small timelike (negative) Q 2 is measured through the 
Dalitz decay it — > e + e~7. In our notations, the usual 
representation of the results is 



K(0,Q 2 



(40) 



where is the experimental pion mass. Then the 
Q 2 -slope given by Eq. ([39]) corresponds to » 0.031. 
This number is not very far from the central values of 
two last experiments, o T = 0.026 ± 0.024 ± 0.0048 [z3, 
= 0.025 ± 0.014 ± 0.026 [75], but the experimental 
errors are rather large. An earlier experiment [76| pro- 
duced = —0.11 ± 0.03 ± 0.08, a result whose central 
value has opposite sign and much larger absolute mag- 
nitude. In the spacelike region, the data are available 
only for the values Q 2 > 0.5 GeV 2 (CELLO [3) and 
Q 2 > 1.5 GeV 2 (CLEO which cannot be treated as 
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very small. The CELLO collaboration [14( gives the value 
Ott = 0.0326 ± 0.0026 that is very close to our result. To 
settle the uncertainty of the timelike data (and also on 
its own grounds), it would be interesting to have data 
on the slope from the spacelike region of very small Q 2 , 
which may be obtained by modification of the PRIMEX 
experiment ffl\ at JLab. 



vanishes for 



B. Large virtualities 

Since ^7{Q,Zo) exponentially ~ e _( ^ 
large Q, we can neglect the first term of Eq. (|3€3[) in 
the asymptotic Q 2 — > oo region. Representing d z ip(z) 
in terms of the &(z) wave function, we write K(Q 2 , Q 2 ) 
as 



K(QlQl) 



so 
2 



J(Qi,z)J{Q 2 ,z)$(z)zdz . (41) 



Our goal is to compare the predictions based on this 
formula with the leading-order perturbative QCD re- 
sults. Note that the situation is different from that with 
the charged pion form factor, where the leading-power 
1/Q 2 pQCD result FpQ cd (Q 2 ) - 2 (o./tt) s /Q 2 [n| 
corresponding to hard contribution is the 0(a s ) cor- 
rection to the soft contribution, for which AdS/QCD 
gives F^ dS/QCD (Q 2 ) -> s /Q 2 |6l|. In that situation, it 
makes no sense to discuss whether pQCD and AdS/QCD 
asymptotic predictions agree with each other numerically 
or not. In general, our AdS/QCD model contains no 
information about hard gluon exchanges of pQCD that 
produce the a s factors. However, the pQCD expression 
for the 7*7* — > 7r° form factor has zero order in a s , so 
now it makes sense to compare the leading-order pQCD 
predictions for this particular (in fact, exceptional) form 
factor with AdS/QCD calculations. 

It is instructive to consider first two simple kinematic 
situations (Q\ = 0, Q\ — Q 2 and Q\ = Q 2 — Q 2 ), and 
then analyze the general case. 



1. One real photon 

Form factor in the kinematics when the virtuality of 
one of the photons can be treated as zero Q 2 « 0, while 
another Q 2 = Q 2 is large was studied experimentally by 
CELLO Q and CLEO [H collaborations. 

In perturbative QCD, the K(0, Q 2 ) form factor at large 
Q 2 is obtained from the factorization formula 



K{0,Q 2 ) = T(Q 2 ,x)<pir(x)dx , (42) 
Jo 

where T{Q 2 ,x) is the amplitude of the hard subprocess 
77* — > qq. The latter, modulo logarithms of Q 2 , has the 
1/Q 2 behavior. In the lowest order, when 

KMa ^'wf/-^ dxs w ,v ' (43) 



the purely 1/Q 2 outcome reflects the large-Q 2 behav- 
ior of the hard quark propagator connecting the photon 
vertices. A particular form of the pion distribution am- 
plitude (DA) is irrelevant to the power of the large-Q 2 
behavior of F(0,Q 2 ), as far as it provides a convergent 
x- integral in Eq. (|43p . The latter requirement is fulfilled, 
e.g., if the pion distribution amplitude ip v (x) vanishes 
at the end-points (x — or 1) as a positive power of 
x(l — x). Whether it vanishes at x = as x, x 2 or ^fx 
does not matter - this would not affect the 1/Q 2 large-Q 2 
behavior of the 77*71° form factor in the lowest pQCD 
order. However, the shape of the pion distribution ampli- 
tude (fir (x) determines the value of the coefficient I v that 
provides the normalization of the C(l/Q 2 ) term. For the 
asymptotic shape <f^ s (x) — 6x(l — x), we have I v — 3 
and ifPQ CD ( as )(0,Q 2 ) = s /Q 2 . 

Take now our extended AdS/QCD model for the 
K(0, Q 2 ) form factor. It gives 



K(0,Q 2 )^^j o 



J{Q lZ )<5>(z)zdz . 



(44) 



At first sight, this expression, though completely different 
analytically from the pQCD formula (|4"2")h has a general 
structure similar to it: the Q 2 -dependence is accumulated 
in the universal current factor 3(Q, z), while the bound 
state dynamics is described by the Q 2 -independent wave 
function $(2). The obvious difference is that the bulk-to- 
boundary propagator J(Q, z) does not have a power be- 
havior at large Q 2 : it behaves in that region like e~® z , co- 
inciding in this limit with zQK\(zQ) = /C(Qz), the free- 
field version of the nonnormalizable mode. The power 
behavior in Q 2 appears only after integration of the ex- 
ponentially decreasing function over z. As a result, only 
small values of z are important in the relevant integral, 
and the outcome is determined by the small-z behav- 
ior of the wave function $(z). As far as Q(z) tends to a 
nonzero value ^(O) when z — > 0, the outcome is the 1/Q 2 
behavior: 



K(0,Q 2 



$(0)SQ 

2Q 2 



d XX 2 K lix) ^^. (45) 



Just like in the case of the (charged) pion EM form fact- 
tor [6l|, the large-Q 2 behavior of K(0, Q 2 ) is determined 
by the value of the $(z) wave function at the origin. 
Note, that this value is fixed: $(0) = 1, which gives 
K(0,Q 2 ) = sq/Q 2 , the result that coincides with the 
leading-order prediction of pQCD for I v = 3, the value 
that is obtained, e.g., if one takes the asymptotic pion 
DA ip^{x) — 6x(l — x). 

Experimentally, the leading-order pQCD prediction 
with I v = 3 is somewhat above the data. However, the 
next-to-leading 0(a s ) pQCD correction is negative, and 
decreases the result by about 15%, producing a satisfac- 
tory agreement. More elaborate fits 22j favor DAs that 
differ from the asymptotic one by higher Gegenbauer har- 
monics x{\— x)C n (2x— 1) with n = 2 and n = 4. Still, for 
all the DAs obtained from these fits, the magnitude of the 
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integral I v is very close to the value I as = 3 given by the does not depend on the shape of the pion DA. In our 



asymptotic DA (see Ref. [22| for details and references). 
Thus, the result of our calculation is in full agreement 
with the magnitude of the leading-order pQCD part of 
the NLO fits of existing experimental data. 



2. Equal virtualities 

Another interesting kinematics is when the photons 
have equal large virtualities, Q\ = Q\ — Q 2 . In this 
case, the leading-order pQCD prediction 



extended AdS/QCD model, we obtain 



K{Q\Q 2 ) 



so 



[J(Q,z)}H(z)zdz . (47) 



Asymptotically, we have 
$(0).s 



K{Q Z ,Q 2 ) 



^ pQCD (Q 2 ,Q 2 ) 



so 



Vir(aj) dx 



Q 2 



dxx 3 [Ki{ X )f 



so 

3Q 2 



(48) 



so 



3 Jo xQ 2 + (1 - x)Q 2 3Q 2 

(46) 



which is the same result as in the leading-order pQCD. 



J 



3. General case 

Finally, let us consider the general kinematics, when Q\ = (1 +uj)Q 2 , and Q\ = (1 — u)Q 2 . The leading-order 
pQCD formula gives in this case 



i fPQC D((1+w)Q 2 )(l _ w)Q 2 ) = ^_^ 1 _ 



(fiir (x) dx 



+ lo(2x-1) 3Q 
while Eq. (|4ip of our AdS/QCD model reduces, for large Q 2 , to 

K({i + lo)q 2 , (i - lo)q 2 ) -> r d x x 3 K^xVT+^^ixVT^) = 4^ /AdS H 



2Q 2 



3Q 2 



(49) 



(50) 



with the function I AdS (cu) given by 



2cj- (1 - lu 2 ) In 



1 -u> 

1 + LO 



(51) 



It is straightforward to check that I AdS (to) coincides with 
the pQCD function calculated for the asymptotic 

distribution amplitude ip as (x) = 6x(l — x). Indeed, using 
the representation 



xKiix) 



we can easily integrate over x m Eq. (|50p to get 



so f°° f°° u 1 u 2 e- Ul - U2 du 1 du 2 
Q 2 Jo Jo u 2 {l + lo) + Ul {\ - uj) 



(52) 



.(53) 



Changing variables u 2 — x\, u± = (1 — x)X and integrat- 
ing over A, we obtain 



K(QlQ 2 2 ) 



•so 



1 6 x(l — x) dx 
3Q 2 Jo 1 +uj(2x- 1) ' 



which coincides with the pQCD formula 
<P-k(x) = 6 x(l — x). 



(54) 
if we take 



Note that the absolute normalization of the asymptotic 
behavior of K(Q\ : Q|) f° r large Qi in our model is fixed 
by the choice k = 2 that allows to conform to the value 
K(Q 2 ,Q 2 ) — 1 corresponding to the QCD anomaly re- 
sult. As we have seen, this choice exactly reproduces also 
the leading-order pQCD result for the equal-virtualities 
form factor K(Q 2 ,Q 2 ). The origin of this rather unex- 
pected result needs further studies. Recall also that in 
pQCD the result for K(Q 2 , Q 2 ) is the same for any pion 
distribution amplitude, while the result for the unequal- 
virtualities form factor K((l + uj)Q 2 , (1 —uj)Q 2 ) depends 
on the choice of the pion distribution amplitude. The 
fact that our AdS/QCD model gives the same result as 
the leading-order pQCD calculation performed for the 
asymptotic distribution amplitude, also deserves a fur- 
ther investigation. 



4- From small to large Q 2 

Both K(0,Q 2 ) and K(Q 2 ,Q 2 ) functions are 
equal to 1 at Q 2 = 0. For large Q 2 , the first one 
tends to so/Q 2 and the second one to so/3Q 2 . The 
question is how these functions interpolate between 
the regions of small and large Q 2 . Long ago, Brodsky 
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and Lepage 
interpolation 



proposed a simple monopole (BL) 



K (0, Q 



1 



(55) 



between the Q 2 = value and the large-Q 2 asymptotic 
prediction of perturbative QCD for K(0, Q 2 ). Later, this 
behavior was obtained within the "local quark-hadron 
duality" approach [27|, [zl], in which K(0, Q 2 ) is ob- 
tained by integrating the spectral density p(s, 0, Q 2 ) = 
Q 2 /{s + Q 2 ) 2 of the 3-point function over the "pion dual- 
ity interval" ^ s < s . The curve for Q 2 K(0, Q 2 ) based 
on Eq. (|44f practically coincides with that based on BL 
interpolation (|55|) . see Figs. [2] and [3] For comparison, we 
also show on Fig. [H the monopole fit AT CLEO (0,Q 2 ) = 
1/T1 + Q 2 /Al) (with A w = 776 MeV) of CLEO data 
fla ]. As we mentioned, an accurate fit to CLEO data [22| 
was obtained in the next-to-leading order (NLO) pQCD, 
with the leading order part of the NLO pQCD fit being 
very close to BL-interpolation curve, and hence, to our 
AdS/QCD result as well. 



K(0,Q 2 ) 




Q 2 (GeV 2 ) 



FIG. 3: Form factor K(0,Q 2 ) in AdS/QCD model (solid 
curve, red online) compared to BL interpolation formula 
(dashed curve, blue online). 



K AdS {Q 2 ,Q 2 )/K LD {Q 2 ,Q 2 ) reaches its maximum value 
of 1.08 for Q 2 ~ 0.3 GeV 2 , then slowly decreasing to- 
wards the limiting value of 1. For large Q 2 , the local 
duality model gives the same result 



Q 2 K(0,Q 2 ) (GeV 2 ) 



K LD (Q\Q*) 




Q 2 (GeV 2 ) 



FIG. 2: Function Q 2 K(0,Q 2 ) in AdS/QCD model (solid 
curve, red online) and in local quark hadron duality 
model (coinciding with Brodsky-Lepage interpolation for- 
mula, dashed curve, blue online). The monopole fit of CLEO 
data is shown by dash-dotted curve (black online). 

In the case of K(Q 2 , Q 2 ) function, our model predicts 
the slope 1.92/m 2 w 2.15/s (twice the slope of K(0 1 Q 2 ), 
see Eq. (|3"9l ), while the local duality model gives [27| 



K hD (Q 2 ,Q 2 ) = l-Q 4 
2Q 2 



dx 



= 1 



o [Q 2 + s x(l - x)} 2 
4 tanh-VV(s + 4Q 2 ) 



so + 4Q 2 y/s {s + 4Q 2 ) 3 

a curve which has the slope 2/sq at Q 2 = 



(56) 



so 

3Q 2 



0(1/Q 4 ) 



(58) 



as our present model (g5]) and pQCD PS]). As a conse- 
quence, our present model produces a curve that is very 
close to the curve based on the local duality model, sec 
Fig.rj 



K{Q\Q 2 ) 




Q 2 (GeV 2 ) 



FIG. 4: Form factor K(Q 2 ,Q 2 ) in AdS/QCD model (solid 
curve, red online) compared to the local quark-hadron duality 
model prediction (dashed curve, blue online). 



V. BOUND-STATE DECOMPOSITION 



^ LD (Q 2 ,Q 2 ) = l-2— +0(Q 4 ) . (57) 
so 

However, higher terms of Q 2 expansion become impor- 
tant for Q 2 as small as 0.01 GeV 2 , where K AdS (Q 2 , Q 2 ) 
becomes larger than K LB (Q 2 ,Q 2 ), and the ratio 



The bulk-to-boundary propagator J7(Q, z) may be 
written as a sum 



J(Q,z)=^ /)2 



^ Q 2 



Ml 



(59) 
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over all vector bound states. For this reason, the form 
factor K(Q\,Q\) also has a generalized vector meson 
dominance (GVMD) representation 



A 7} ._k + Br, 



K(QlQl) -gg (1 + Qf/ ^ )(1 + Q 



(60) 



where A n ,k$ come from the first (surface) term in 
Eq.®, 



.4 



n.k 



4*(z ) 



(61) 



while 's are obtained from the second term and are 
given by the convolutions 



B„ 



2o70 : n70,fc J\ (lO,n)J\ (lO,k) 



x / V'(z) Jii^nz/zo) J 1 {j ,kz/z )z 2 dz . (62) 



Let us study the structure of the bound-state decompo- 
sition in two most interesting cases: for K(0, Q 2 ) and 
K(Q 2 7 Q 2 ). 



This relation may be directly obtained from the formula 



n~[ 70,n^l 2 (70,n) 



1 



(69) 



The bound state decomposition of K(0, 0) looks like 

K(0, 0) = 1 = 0.9512+0.0408 + 0.0446 - 0.0753 (70) 
+0.0764 - 0.0736 + 0.0703 + . . . . 

There is a strong dominance of the lowest vector state, 
while each of the higher states is suppressed by more than 
factor of 10. The slow convergence of higher terms is due 
to A n terms proportional to $f(zo) m 0.14. For large n, 
one can approximate A n w ^f(zo)(—l) n y/2/n. 

Integrating by parts in Eq. (|66|) gives a representation 
directly for the total coefficient 

2 f z ° 

A n + B n = 9 t9 , / *(z) Jo(7o,n«Ao) zdz , 

'o 

(71) 



z o J i(7o,-. 



that is related to the expansion of the pion wave function 
^(z) over the </> ra -functions (JSJ of vector meson bound 
states: 



A. One real photon 

To study the bound-state decomposition of K(0, Q 2 ), 
we write the basic expression 

K(0, Q 2 ) = *(z )J(Q, z ) - / J(Q, z) V'(z) dz (63) 



and use GVMD representation ([59} for J(Q,z). This 
gives 



where 



2*(z ) 



70,n-M70,n) 
and BnS coincide with the coefficients 
2 



(64) 



(65) 



B n — — - 



%q To, n JKlO-n) Jo 
for the expansion 



V'(z) J 1 (- f0 ,nz/z )zdz (66) 



1 ^ 

*'(«) = V S„ 7o ,„ Ji (7o,n«M>) (67) 

n— 1 

of the pion wave function \l/'(z) over the ipX(z)/z wave 
functions (T5T) of vector meson bound states. In particular, 



oo oo 



^(0,0) = E^« + E B " 

n— 1 n— 1 

= tf (z„) + [* (o) - = i • (68) 



(A„ + S n ) J ("/ 0: nZ/z ) 



(72) 



Using it, one obtains again that if(0,0) = ^(0) = 1. 
The slow convergence of higher terms in Eq. (|70j) is now 
explained by the necessity to reproduce the finite value 
of ^f(z) at z — zq by functions vanishing at z = zq. 

The slope of K(0, Q 2 ) at Q 2 = is given by the sum of 
(A n + B n )/M 2 coefficients, which converges rather fast: 

K(0,Q 2 ) = l-^(o.9512 + 0.0077 + 0.0034 - 0.0031 



0.0020- 0.0013 + 0.0009 



(73) 



and the contribution of the lowest state completely dom- 
inates the outcome. 

Each term of the GVMD expansion behaves like 
1/Q 2 at large Q 2 . In particular, the lowest-state con- 
tribution behaves like 0.95 m 2 p /Q 2 « (0.57 GeV 2 )/Q 2 . 
We also obtained above that K(0, Q 2 ) behaves asymp- 
totically like s Q /Q 2 « (0.67GeV 2 )/Q 2 . The two scales 
are not so different, and one may be tempted to spec- 
ulate that the large-Q 2 behavior of K(0, Q 2 ) also re- 
flects the dominance of the lowest resonance. However, 
the coefficient of 1/Q 2 is formally given by the sum of 
(A n + B n )M 2 terms, which does not show good conver- 
gence even after 7 terms are taken: 

OO s 

^T(A n + B n )M 2 = m 2 J 0.951 + 0.215 + 0.577 

n=l *- 



- 1.811 + 2.945 - 4.158 + 5.473 



(74) 
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A simple comment is in order now: within the 
AdS/QCD model [43[ the p-meson mass is determined by 
the "confinement radius" zq, while the scale so = 8n 2 f 2 
is basically determined by the chiral symmetry break- 
ing parameter a (see Eq. (II 6[) and preceding discussion) . 
Calculationally, the coefficient so of 1/Q 2 asymptotic be- 
havior was determined solely by the magnitude of the 
pion wave function ^'{z)/z at the origin. Furthermore, 
it was legitimate to take the free-field form of the vec- 
tor bulk-to-boundary propagator in our calculation, i.e., 
no information about vector channel mass scales was in- 
volved. 

Moreover, one may write the bound-state decompo- 
sition for the J(Q, zq) function. Again, each term of 
such a decomposition has 1/Q 2 asymptotic behavior, 
while J(Q, zq) exponentially decreases for large Q. In 
fact, the formal sum J2 n A n M 2 in this case diverges like 

£ n (-l)«„3/2. 

Summarizing, the 1/Q 2 asymptotic behavior of 
K(0, Q 2 ) has nothing to do with the fact that the contri- 
bution of each particular bound state has 1/Q 2 behavior. 
If, instead of $(z), one would take a function with ~ z A 
behavior for small z, one would still be able to write the 
GVMD representation for such a version of K(0, Q 2 ), but 
its asymptotic behavior will be ~ 1/Q 2+A . 



B. Two deeply virtual photons 



K(Q 2 ,Q 2 ) behaves like 1/Q 2 for large Q 2 . This result 
was a consequence of two features of the form factor inte- 
gral pi)) . The first is the fact that the bulk-to-boundary 
propagator J{Q, z) behaves like e~® z for large Q. This 
is a very general property: in this limit J(Q, z) should 
coincide with its free-field version IC(Qz) — zQKi(zQ). 
The second feature is that the pion wave function is 
finite at the origin, which follows from the basic formula 
1(15)1 that defines / w . 

Hence, to qualitatively understand the mechanism that 
produces 1/Q 2 result from the doubly- infinite sum of 
1/Q 4 terms, one may consider a simpler "toy" model that 
also has these general properties, but allows to analyti- 
cally calculate integrals that determine the coefficients in 
Eq. ((60)) . In particular, an explicit result for K(Q 2 ,Q 2 ) 
may be obtained if we take the soft-model expression [691 ] 



J s (Q,z) = a 



exp 



2 2 
■ K Z 



dx (75) 



for the bulk-to-boundary propagators (with k being the 
oscillator parameter and o = Q 2 /4k 2 ) and 

y s (z) = e - K2 * 2 (76) 

for the pion wave function. This model has the required 
properties, namely, J s (Q,z) approaches the free-field 
function K.{Qz) for large Q 2 , while $(0) is finite. 
Calculating the integral 



Each term of the bound-state decomposition ((60)) for 
K(Q 2 , Q 2 ) has 1/Q 4 behavior. Thus, in the case of strong 
dominance of a few lowest states, one would expect 1/Q 4 
large-Q 2 behavior of this function. 

However, as we already obtained, the function 



K S (Q(, Qi) = 2k z / J s (Qi,z) J S (Q 2 , z) e~ K z zdz 



gives 



(77) 



K s (Ql Ql) = E i z )( a \ rrw 7 u \ i^n > ( ?8 ) 

(oi + n){ai +n + l) (a 2 + n)(a 2 +n + 1) 

where a ± = Qj/M 2 , a 2 = Qj/M 2 and M = 2 K is the mass of the lowest bound state. The spectrum corresponding 
to J S (Q, z) is given by = mM 2 , with m = 1, 2, . . . , and Eq. 1(75)) explicitly displays the bound state poles. For 
large Q 2 ,Q 2 , each term of this sum behaves like l/QiQ2- However, taking ai = a 2 = a ^> 1 gives 

r°° rln 1 M 2 

K%Q 2 ,Q 2 )^a 2 ^— = ±=M. (79) 
Jo (n + a) 4 3a 3Q Z 

Thus, the conversion from the 1/Q 4 asymptotics of individual terms to the 1/Q 2 asymptotics of the sum is due to 
nondecreasing 0(nP) behavior of the coefficients accompanying n th term of the sum. In other words, transitions 
involving higher bound states are important, i.e., the pole decomposition is far from being dominated by a few lowest 
states. 



C. Structure of two-channel pole decomposition 



However, to make specific statements about the transitions, one should realize that Eq. (|78p does not have the 
form of Eq. (|60p . In particular, it is not a double sum, and having a sum over a single parameter implies that the 
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summation parameters n and k in the double sum representation are correlated. A simple inspection of Eq. (|78|l shows 
that either k = n or k = n±l. Furthermore, the representation (|75|) has Q\Q\ factor in the numerator, which should 
be canceled against the denominator factors to get an expression in which Q\ and Q\ appear in denominators only. 
The easiest way to obtain the desired GVMD-type expansion for K S (Q 2 , Q 2 ) is to use another representation [69[ 



J s (Q,z) = k 2 z 2 I exp 
'o 



1-x 



2 2 
K Z 



x a dx 



(80) 



for the bulk-to-boundary propagators. Then 



K*(QlQt) = J2 



n=0 



(n + l)(n + 2) 



(ai + n+ l)(ai + n + 2)(a 2 + n + l)(a 2 + n + 2) 



(81) 



Now, each term of the sum decreases as 1/QfQf, but the sum may be rewritten as 



1 



2 /M 2 l + Q 2 2 /M 2 +1 1 + Ql/M^J' 



(82) 



i.e., the coefficients B s n . of the bound-state expansion ([SO]) in this model are given by 



(83) 



(there is no need to add surface terms in this model since ^ s (oo) = 0, and hence A s n k = 0). 

I ' 



When Q\ = 0, only the n = 1 term contributes, and 
K(Q 2 , 0) in this model is formally dominated by the low- 
est resonance: 



1 



1 + Q 2 /M( 



(84) 



k\n 


1 


2 


3 


4 


5 


1 


0.7905 


0.0272 


-0.1331 


0.0562 


-0.0275 


2 


0.0272 


0.2484 


-0.0423 


-0.0608 


0.0287 


3 


-0.1331 


-0.0423 


0.1624 


-0.0367 


-0.0403 


4 


0.0562 


-0.0608 


-0.0367 


0.1199 


-0.0303 


5 


-0.0275 


0.0287 


-0.0403 


-0.0303 


0.0951 



TABLE I: Coefficients B n t in the hard-wall model. 



This fact may create an impression that the ir° — * 77 
decay in this model is dominated by the pu> intermediate 
state. However, the outcome that the sums of bracketed 
terms are zero for n > 2 is a result of cancellation of the 
contribution of a diagonal transition that gives 2, and 
two off-diagonal transitions, each of which gives —1. 

In fact, the coefficients B s n n — 2 of diagonal transitions 
do not depend on n, and their total contribution into 
K S (Q 2 , Q 2 ) diverges. On the other hand, the coefficients 
B^ n+1 = —1 and B^ +l n = —1 of subdiagonal transi- 
tions are negative, and also do not depend on 



total contribution into K S {Q\ 1 Q\) of each of k = n 



The 
1 

or k = n — 1 off-diagonal transitions also diverges. In 
such a situation, claiming the dominance of the lowest 
states contribution makes no sense. 



In the hard-wall model, the diagonal coefficients B n ^ n 
decrease with n, and they are visibly larger than the 
neighboring off-diagonal ones (see Table I). Thus, one 
may say that, for small Ql,Q 2 , the value of K{Q\, Q|) 
is dominated by the lowest bound states (the coefficients 
A n n also decrease with n, see Table II, asymptotically 
they behave like 2^{zo)/n), However, the large-Q^Q^ 



k\n 


1 


2 


3 


4 


5 


1 


0.3641 


-0.2420 


0.1935 


-0.1658 


0.1474 


2 


-0.2420 


0.1609 


-0.1286 


0.1102 


-0.0980 


3 


0.1935 


-0.1286 


0.1028 


-0.0881 


0.0783 


4 


-0.1658 


0.1102 


-0.0881 


0.0755 


-0.0671 


5 


0.1474 


-0.0980 


0.0783 


-0.0671 


0.0597 



TABLE II: Coefficients A n k in the hard- wall model 



expansion is given by 
K(n 2 n 2, Mt ™ M 2 M 2 



(85) 



i.e., one should deal with the coefficients 



B,,^S=B^^ (86) 



Mf 



7o 4 ,i 



Mj 4 



the lowest of which are given in Tables III and IV. Again, 



k\n 


1 


2 


3 


4 


5 


1 


0.7905 


0.1433 


-1.7236 


1.3514 


-1.0593 


2 


0.1433 


6.8952 


-2.8875 


-7.6956 


5.8234 


3 


-1.7236 


-2.8875 


27.2323 


-11.4306 


-20.1136 


4 


1.3514 


-7.6956 


-11.4306 


69.3299 


-28.0684 


5 


-1.0593 


5.8234 


-20.1136 


-28.0684 


141.2499 



TABLE III: Coefficients B nyk M%Ml/Mf in the hard- wall 
model. 



k\n 


1 


2 


3 


4 


5 


1 


0.3641 


-1.2751 


2.5056 


-3.9868 


5.6816 


2 


-1.2751 


4.4655 


-8.775 


13.9624 


-19.8979 


3 


2.5056 


-8.775 


17.2437 


-27.4374 


39.1011 


4 


-3.9868 


13.9624 


-27.4374 


43.6572 


-62.216 


5 


5.6816 


-19.8979 


39.1011 


-62.216 


88.6642 



TABLE IV: Coefficients A„, k M%Mi/Mf in the hard-wall 
model 



the coefficients increase with n and k producing diver- 
gent series, just like in the toy soft-like model. Note, 
that asymptotically the sum of A-type terms gives a con- 
tribution exponentially decreasing with Q\ and/or Q2, 
i.e., much faster than the ~ 1/Q 4 asymptotic behavior of 
each particular transition. On the other hand, the sum 
of .B-type terms gives a contribution that has ~ 1/Q 2 
asymptotic behavior, i.e., it drops slower than 1/Q 4 . 

The convergence situation may be different in the real- 
world QCD, in which higher resonances are broad, with 
the width increasing with n (or k). Then the diagonal 
and neighboring non-diagonal transitions strongly over- 
lap for large n and may essentially cancel each other. 



VI. SUMMARY 

At the end of the pioneering paper [43j , it was indicated 
that one of the future developments of the holographic 
models would be an incorporation of the 5D Chern- 
Simons term to reproduce the chiral anomaly of QCD. 
However, only relatively recently Ref. [HtJ discussed a 
holographic model of QCD that includes Chern-Simons 
term (see also [79|) and, furthermore, extends the gauged 
SU(2) L ® SU(2) R flavor group to U(2) L <g> U(2) R . 

In the present paper, we develop an extension of the 
AdS/QCD model, similar in form to that proposed in 
[67| . but adjusted to study the anomalous coupling of 
the neutral pion to two (in general, virtual) photons. 
The additional part of the gauge held in the 5D bulk 
is associated with the isoscalar vector current (related 
to w-like mesons). The Chern-Simons term allows to 
reproduce the tensor structure of the anomalous form 
factor F 1 , J » n o(Ql, Q\). To exactly reproduce the QCD 
anomaly result for real photons, we added contributions 
localized at the IR boundary z = z$, and then studied 
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the momentum dependence of the _F! y » 7 » 7r o (Q 2 , Q 2 .) form 
factor in our model. 

In particular, we calculated the slope of the form 
factor with one real and one slightly off-shell photon. 
Our result f» 0.031 for the parameter of the usual 
F r? ,^ (0, Q 2 ) = F^.^o (0, 0)(1 - a7T Q 2 /ml) low-Q 2 ex- 
perimental representation of the data is very close to the 
value a n = 0.0326 ± 0.0026 obtained by CELLO collabo- 
ration [LJ] from spacelike Q 2 measurements, and rather 
close to the central values ~ 0.024 of two most recent 
experiments |74l |75| for timelike Q 2 . 

Although the holographic model is expected to work 
for low energies, where QCD is in the strong coupling 
regime, we found it interesting to investigate the behav- 
ior of the model form factor also in the regions where at 
least one of the photon virtualities is large. For the case 
with one real and one highly virtual photon, we demon- 
strated that our AdS/QCD result is in full agreement 
with the magnitude of the leading-order part of the next- 
to-leading-order pQCD fits of existing experimental data. 
In the kinematics when both photons have equal and 
large virtualities we obtained the same result as in the 
leading-order pQCD. Finally, we considered the general 
case of unequal and large photon virtualities. In this case, 
the form factor has a nontrivial dependence on the ratio 
Q1/Q2 °f photon virtualities. Our calculation shows that 
the hnal result of our AdS/QCD model analytically co- 
incides with the pQCD expression calculated using the 
asymptotic distribution amplitude ip^{x) = 6x(l — x). 
This result is rather unexpected, because initial expres- 
sions for the form factor have very different structure. It 
should be noted that the absolute normalization of the 
form factor K{Q\ : Q\) in our model is fixed by adjust- 
ing its value to K(0, 0) = 1 at the real photon point, 
which allows to conform to the QCD axial anomaly. The 
outcome that this choice exactly reproduces the leading- 
order pQCD result for the equal-virtualities form factor 
K(Q 2 , Q 2 ) needs further studies, as well as our result that 
the w-dependence of the unequal-virtualities form factor 
K ((1 +ui)Q 2 , (1 — lo)Q 2 ) coincides with the leading-order 
pQCD result derived by assuming the asymptotic shape 
for the pion distribution amplitude. 

The bulk-to-boundary propagators entering into 
AdS/QCD formulas for form factors have a generalized 
vector-meson-dominance (GVMD) decomposition. As a 
result, the form factors also can be written in GVMD 
form. We studied the interplay between the GVMD 
decomposition of form factors and their behavior for 
large photon virtualities. In the case of one real pho- 
ton, the function K(0,Q 2 ) asymptotically behaves like 
1/Q 2 . However, we demonstrated that this behavior has 
nothing to do with the fact that each term of the GVMD 
expansion for K(0,Q 2 ) also behaves like 1/Q 2 for large 
Q 2 . In fact, a formal GVMD expression for the coeffi- 
cient of the 1/Q 2 term diverges. When both photons are 
highly virtual, each term of the GVMD expansion for 
K(Q 2 ,Q 2 ) behaves like 1/Q 4 , while K(Q 2 ,Q 2 ) has 1/Q 2 
asymptotic behavior. Thus, we observe that only in the 
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region of small photon virtualities it makes sense to talk 
about dominating role of the lowest states. In partic- 
ular, for real photons, when Q\ = Q\ = 0, the lowest 
("/9cl>7t") transition amplitude contributes 1.15 into the 
K(Q, 0) = 1 value, the excess being primarily cancelled 
by the neighboring non-diagonal transitions. 
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